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It is shown that the voltage dependence of the tunneling conductance between two lightly doped 
semiconductors, which are separated by an large area tunneling barrier, can reveal the high energy 
part of the Coulomb gap if the barrier is thick enough. At the barrier thickness smaller than 
average distance between impurities no Coulomb gap feature can be found. This happens because 
such tunneling is sensitive to very rare shortest pairs of occupied and empty states localized at 
opposite sides of the barrier, whose density of states in this limit has no Coulomb gap. Small area 
tunneling contacts are also discussed. It is shown that the tunneling conductance of a point-like 
contact exponentially grows with the applied voltage. This dependence does not permit a direct 
measurement of the Coulomb gap. 



The density of states of one-electron (or more generally 
charged) excitations, g(e), of a three-dimensional sys- 
tem with localized electrons is known to have the 
Coulomb gap near the Fermi level, 



g(e) = (3/7r) K 3 e' 



(1) 



where e is energy calculated from the Fermi level, e is the 
proton charge k is the dielectric constant of the semicon- 
ductor. Eq. (0) is valid while e < A ~ e 2 /kR, where R is 
the average distance between impurities in the semicon- 
ductor, A is the width of the Coulomb gap and simulta- 
neously the width of the peak of the density of states of a 
lightly doped and moderately compensated semiconduc- 
tor. 

The Coulomb gap is known to dramatically change DC 
variable range hopping conductivity leading to the Efros- 
Shklovskii law ||,| 



a = ( T cxp[-(To/T) 1 / 2 ] 



(2) 



Here To = Ce 2 /na, where a is the decay length (effective 
Bohr radius) of impurity states and C is the numerical 
coefficient. 

On the other hand, the phononlcss high frequency hop- 
ping conductivity is not sensitive to the Coulomb gap, 
because it is determined by very close pairs of empty 
and occupied localized states QIJ. 

An appealing method of direct observation of the 
Coulomb gap is low temperature tunneling into a semi- 
conductor from a metallic electrode separated by a tun- 
neling barrier. If the barrier is thick enough, screening 
of the Coulomb interaction in the semiconductor by the 
metal can be neglected. In this case, the density of local- 
ized states in the semiconductor can be measured study- 
ing variation of the tunneling conductance as a function 
of the applied voltage. Several works [^|-(6) claim obser- 
vation of the Coulomb gap in such a way, although, it is 
not totally clear why in this cases screening by the metal 
was weak enough. 



Recently, tunneling through the break junction in a 
germanium crystal with localized states was studied M. 
The tunnel conductance showed a minimum at zero bias 
voltage which was interpreted as the Coulomb gap. 

In this paper, we consider theoretically tunneling be- 
tween two semiconductors, where states are localized. 
Both of them can be three-dimensional, so that a bar- 
rier separates two half-spaces. In the two-dimensional 
case, when electrons occupy only one plane, the barrier 
can be a stripe separating two half-planes filled by lo- 
calized two-dimensional electrons. Another system, with 
the same physics consists of two parallel layers of local- 
ized two-dimensional electrons separated by a thin bar- 
rier layer. For certainty, we concentrate here on the 
three-dimensional case and assume that parameters of 
semiconductors on both sides of the barrier are exactly 
the same. 

Tunneling between semiconductors with discrete local- 
ized states has to be due to phonon assisted processes, be- 
cause probability to find exactly resonant pairs of states 
on both sides of the barrier is infinitesimally small. In 
this paper we concentrate on low temperature tunneling 
which goes with phonon emission only. 

The main goal of this paper is to understand how the 
current-voltage characteristics of such a tunneling con- 
tact depends on the barrier thickness d. If the barrier is 
thin, d <C R, the distribution of electrons in the ground 
state and, therefore, the Coulomb gap of the density of 
states is not affected by the barrier. We assume that 
density of states in both semiconductors is still given by 
Eq. (gj) even if d > R @. 

Below we assume that the tunneling barrier is high 
enough so that all resistance is determined by hops from 
occupied states on the left side of the barrier to empty 
states on the right side. In other words, if we introduce 
the Miller-Abrahams resistance network || connecting 
all the localized states inside each semiconductor and be- 
tween them we immediately get an enormous difference 
between these two types of resistors. Resistors crossing 
the barrier have much larger resistance and, therefore, 
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they determine the tunneling conductance between two 
semiconductors. One can imagine that the bulk of both 
semiconductors is a superconductor and the tunneling 
conductance is due to parallel resistors which cross the 
barrier. Then it is clear that at large enough d, the num- 
ber of these parallel resistors and, therefore, the tunnel- 
ing current, /, is proportional to the product of number 
of localized states with energy from < e < ell/ 2 on the 
left side and number of empty states in the interval rom 
> e > ~eU/2. Each of these numbers is proportional 
to U 3 , so that the tunneling current is equal 

I = A{U){eUf. (3) 

where coefficient A(U) depends on U through the phonon 
energy and is specified below. We show below that at a 
finite barrier thickness d the current-voltage character- 
istic is more complicated than Eq. ([}]), namely, that at 
U + e 2 1 nd < A the current is 

I = A{U)eU{eU + e 2 / ndf . (4) 

In order to derive Eq. ([!]), let us following Ref. || in- 
troduce such a function F(fl, r) that F(fl, r)4:TTr 2 drdQ is 
the concentration of pairs of occupied and empty states 
with the distance between them (the arm of the pair) in 
the range between r and r + dr and with the energy re- 
quired for electron transfer in the interval between SI and 
SI + dfl. This function was calculated in Ref. |3) in the 
following way: 

p oo f oo 2 

F(n, r) = / / de 1 de 2 g{e 1 )g{e 2 )6{e 1 - £2 - — - ft) 
Jo Jo Kr 

= B K V 12 (0 + e 2 / K r) 5 , (5) 

where B = 3/(107r 2 ). Using this function one can eas- 
ily calculate the tunneling current /. Indeed all pairs 
crossing tunneling barrier can contribute to /. They 
all work in parallel because as we said above the cur- 
rent reaches their ends almost without resistance. Let us 
write the rate of phonon assisted tunneling hops across 
the barrier with emission of a phonon with energy SI as 
v(Q) = vq(Q) exp(— 2d /ao). At low enough SI only acous- 
tic phonons are involved. For emission acoustic phonons 
via deformational potential interaction vq(Q) oc SI, when 
SI < hs/a, where s is the sound velocity. At larger SI the 
rate i/(Sl) sharply decreases with SI because wavelength 
of the phonon becomes smaller than a (see, for example, 
Ref. §). 

Now we can identify optimal pairs which determine 
the contact conductance. They all have arms which are 
very close to r = d. Indeed, the barrier part of the arm 
of such a pair can not be longer than d + oq, while the 
semiconductor part of the arm can not be longer than a. 
Otherwise, a pair has exponentially smaller than optimal 
ones conductance and should be disregarded. Thus, the 
occupied state of an optimal pair should be in the layer 



with thickness a near the plane of the left semiconductor 
barrier junction. If area of the junction is S the volume 
available for the occupied pair is Sa. For a given oc- 
cupied state the empty state of an optimal pair should 
reside in the cylinder with the height a and with the base 
with the shape of the disk of the radius (dao) 1 ^ 2 drawn 
on the plane of the right barrier-semiconductor junction 
around the center which exactly opposes the occupied 
state on the left junction. The transfer energy, SI, should 
be in the interval < SI < ell. Otherwise, the pair does 
not consist of occupied and empty states and does not 
contribute to the tunneling current. Thus, total num- 
ber of optimal pairs is N = F(eU, d)eUSa 2 aod. For the 
tunneling current Nev we arrive at Eq. (||) where 

A(U) = BK 6 e- n (Sa 2 aod)u (eU) exp(-2d/a ). (6) 

Looking at Eqs. (|), (§), we see that the current- 
voltage characteristics has two voltage dependent factors. 
One of them, the rate of emission of phonons, fo(eU) is 
pretty well known. The ratio I(U)/Uv(U) gives the en- 
ergy dependence apparent join density of states of the two 
semiconductors. At small U this joint density of states 
does not depend on U. If d 3> R, only at large enough 
voltages, ell > e 2 /nd, one can recover the high energy 
part of the joint density proportional to U 5 and directly 
indicating existence of the Coulomb gap. On the other 
hand, if d < R, the joint density of states is constant at 
all elf < A. At larger U the tunneling current decreases 
due the decrease of the density of states of the classical 
impurity band at energies larger than A Q. Thus, at 
d < R one can not see the Coulomb gap at all. 

The reason for the absence of [/-dependence of the joint 
density of states at small d is in the strong Coulomb in- 
teraction of the electron transferred to the right side of 
the barrier with the hole remaining on the left side. The 
same effect leads to the enhancement of the phononless 
high frequency conductivity || . Both phenomena are de- 
termined by rare, isolated from each other pairs of states. 
When comparing Eq. (Q) with formulas for the high fre- 
quency conductivity one can easily find straightforward 
analogies if eU and d of this paper are replaced by ftio and 
by the arm, r u , of pairs which absorb radiation quanta 
hio. 

Our theory is applicable only to a barrier which height 
and thickness are uniform along a large area, S, of the 
contact which has many optimal pairs. If area of the 
barrier, S, is so small that in average there is less than 
one optimal pair, the situation becomes mesoscopic. This 
happens at S ~ S c , where 

S c = K-*e a {eU)-\eU + e 2 / nd)- 5 a" 2 (aod)- 1 . (7) 

At S <C Sc the tunneling conductance is determined by 
the best pair which can be found in this area of the bar- 
rier. Such a pair is longer than d by a random length 
Sr. Therefore, tunneling conductance of a barrier with 
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S <C S c has an additional small factor exp(— 2Sr/a). Of 
course, in such conditions the conductance of the tunnel- 
ing contact strongly fluctuates from a sample to sample. 
This is why we call this case mesoscopic. 

Let us consider how growing voltage affects the conduc- 
tance of a tunneling contact with a small area S. With 
growing U new pairs with smaller Sr become available 
and the exponentially small factor, exp(— 2Sr/a), grows. 
Eventually it reaches unity when Sr becomes as small as 
a. Simultaneously, as we see from Eq. (Q) , critical area S c 
decreases with growing U and eventually approaches S. 
At this point we return to the conductance of a macro- 
scopic barrier. This means that while S c ^ S the tunnel- 
ing conductance grows exponentially with voltage U sim- 
ulating an the joint density of states and the one-electron 
density of states g^e) exponentially growing with energy. 
This phenomenon is, of course, different from a direct 
measurement of the Coulomb gap. Thus, while for the 
uniform barrier the Coulomb gap can be revealed at large 
enough d, the point contact does not directly reveal the 
Coulomb gap. 
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